Acoustic multipole sources from the Boltzmann equation 
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By adding a particle source term in the Boltzmann equation of kinetic theory, it is possible to 
represent particles appearing and disappearing throughout the fluid with a specified distribution of 
particle velocities. By deriving the wave equation from this modified Boltzmann equation via the 
conservation equations of fluid mechanics, multipole source terms in the wave equation are found. 
These multipole source terms are given by the particle source term in the Boltzmann equation. To 
the Euler level in the momentum equation, a monopole and a dipole source term appear in the wave 
equation. To the Navier-Stokes level, a quadrupole term with negligible magnitude also appears. 



I. INTRODUCTION 



Acoustic multipoles are oscillating sources that emit 
acoustic fields of different directivities. These sources can 
be either point sources, localized at single points in space, 
or they can be distributions throughout the medium. The 
the first three orders of multipoles are the most well- 
known: Monopoles, dipoles, and quadrupoles at zeroth, 
first, and second order, respectively. 

When these three types of multipole sources appear as 
source terms in the wave equation, they usually originate 
from terms in the conservation equations of fluid mechan- 
ics. For instance, monopoles are linked to a source term 
in the mass conservation equation (also known as the con- 
tinuity equation), which represents mass appearing and 
disappearing throughout the fluid as a function of time. 
This can model pulsations of small bodies throughout the 
fluid Q. 

However, an alternative approach is to add a particle 
source term in the Boltzmann equation, which is more 
fundamental than the fluid conservation equations which 
can be derived from it. Adding such a source term to the 
Boltzmann equation allows specifying the velocity distri- 
bution of particles that appear and disappear through- 
out the fluid. This approach is therefore similar to, but 
more general than, the aforementioned method of adding 
a mass source term, and can therefore model more gen- 
eral vibrations of small bodies in the fluid. 

Such a particle source term was recently examined for 
the lattice Boltzmann method Q, a computational fluid 
dynamics method based on the fully discretised Boltz- 
mann equation. It was found that such an approach re- 
sults in a wave equation with non-vanishing monopole, 
dipole, and quadrupole source terms. This article will 
similarly examine particle source terms, but in the clas- 
sic non-discretised Boltzmann equation. 



II. ACOUSTIC MULTIPOLE SOURCES 

Mathematically, multipoles are related to source terms 
in the wave equation, 



l__d 



- V' z ) P (x, t) = T (x, t) + ^-Ti(x, t) 



if- 



dxidx 



■T i3 {x,t) + . . 



(1) 



Here, p is the pressure and cq is the speed of sound. The 
terms on the right-hand side are multipole source terms. 

This article makes use of the index notation commonly 
used in the field of fluid mechanics. In this notation, a 
single index indicates a generic vector element (e.g. T, 
could be T x , T y , or T z ), and multiple indices (as in Ty or 
aibj) indicate generic tensor elements. Repeating indices 
within a single term implies summation over all possible 
values of that index. For example, a^fe; = a x b x + a y b y + 
a z b z = a b, and dTi/dxi = 8T x /dx+dT y /dy+dT z /dz = 
V • T. When indices repeat in this way, the letter used 
is arbitrary, so that a^fr; = akbk- 

The general three-dimensional solution to ((T]) is given 
by an integral over the entire volume of the source terms 
on the right-hand side [l|, 
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Thus, T (x,t) indicates the instantaneous monopole 
strength, Ti(x,t) the z-dipole strength, and the ij- 
quadrupole strength. 

As mentioned above, monopole sources can be modeled 
by adding a mass source term to the mass conservation 
equation, 
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p being the mass density, u the fluid velocity, and Q(x, t) 
the instantaneous mass flux. Dipoles typically originate 
from the force term in the momentum conservation equa- 
tion. To the Euler level, this equation is 
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where F represents body forces. Finally, quadrupoles 
typically originate from the nonlinear term in (j?)). 

The multipole terms in ((1) are usually related to the 
terms in the conservation equations as |l| 
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III. THE BOLTZMANN EQUATION 

The Boltzmann equation describes motion of a gas at 
a finer level of detail than the fluid conservation equa- 
tions. In this discussion we shall restrict ourselves to 
its very basics. More details can be found in the litera- 
ture [3|-l5j. The equation evolves the particle distribution 
function, f(x,£,t), which may be seen as a double den- 
sity in both physical space and particle velocity space. 
Thus, it describes the density of particles with position 
x and velocity £ at time t. 

The familiar macroscopic quantities can be recovered 
as moments of /, i.e. by weighting with some function 
and integrating over the entire velocity space. The mass 
density and momentum density are found as the zeroth- 
and first-order moments, 



p(x,t) = J f(x,$,t) d& 
pu(x,t) = J Zf(x,Z,t)d£. 



(5a) 
(5b) 



Neglecting body forces and using the BGK collision 
operator the Boltzmann equation is 
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six, £, t) is the aforementioned particle source term which 
is central to this article. As the left-hand side of the equa- 
tion is a standard advection equation, s(x, £, t) describes 
the rate at which particles are added into the f(x,£,t) 
distribution. The final term is the BGK collision oper- 
ator, which models collisions between particles as a re- 
laxation with a characteristic relaxation time r to the 
equilibrium distribution function, 
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with p and u found from In this article, p is approx- 
imated by ||HJ| . As the quantities of mass and momentum 
are conserved in collisions, substituting /w for / in (J3J) 



must give the same moments. As a result of this, the 
BGK collision operator conserves both mass and momen- 
tum. 

The BGK collision operator is a far simpler model of 
collisions than Boltzmann's original and more accurate 
collision operator. That simplicity comes with draw- 
backs, chiefly that the BGK operator slightly mispredicts 
the Prandtl number Q. This dimensionless number re- 
lates the transport coefficients in the momentum equa- 
tion (viscosity) and the energy equation (conductivity). 
However, this will not matter as we will neglect the effects 
of conductivity in this article. 

From this point on we will assume that the macroscopic 
variables fluctuate only slightly around rest state values 
p = po, p = po, u = 0. This allows us to linearise 
subsequent equations in this article, which is in keeping 
with the usual assumptions in acoustics. 

The pressure p in Q is approximated using the com- 
mon isentropic relation [l|, |7| 
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7 = (d+2)/2 being the adiabatic index determined by the 
degrees of freedom d of the molecules that make up the 
gas Q • In this way, we include the effect of equipartition 
of energy between translational and inner (i.e. rotational 
and vibrational) degrees of freedom, in the limit of rapid 
energy transfer. This relation leads to an ideal speed of 
sound cq given by 
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It will be useful to introduce an abbreviated notation 
for the moments of the particle source term s, 



S (x,t) = J s(x,t,t)d£, 
Si{x,t) = J Zis(x,£,t) d£, 
Sij(x,t)= f Zi£js(x,£,t)d£, 



(10a) 
(10b) 
(10c) 



and so forth. So(x,t) represents the instantaneous mass 
flux of the particles at x, Si is associated with odd sym- 
metries of s in velocity space, and Sij is similarly associ- 
ated with various even symmetries. 



IV. FLUID CONSERVATION EQUATIONS 

It is possible to find the conservation equations of fluid 
mechanics from the Boltzmann equation ^j. To find the 
Euler equations, we could simply take the zeroth and 
first moments of ((6]) under the assumption that / ~ J 1 - '. 
However, to find the momentum equation to the Navier- 
Stokes level, so that it includes the stress tensor term, we 
must resort to the Chapman- Enskog expansion. This is a 
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technique used to derive the fluid conservation equations 
from the Boltzmann equation. It is discussed throughout 
the literature with varying approaches and varying levels 
of complexity [H-d, H> Q • In the following derivation, we 
use a moment- based approach 

Two mathematical techniques are used in this expan- 
sion. First, the distribution function / is approximated 
as a perturbation expansion around equilibrium /(°) in a 
smallness parameter e. Second, a multi-scale expansion 
of time is performed. In mathematical notation, 

+ e\f^ + ... , 
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The smallness parameter e is associated with the dimen- 
sionless Knudsen number Kn = l m f p /L, relating the mean 
free path Z m f p in the gas to a macroscopic length scale L. 
Thus, is of one order higher in the Knudsen num- 

ber than A dimensional analysis Q reveals that 

the relaxation time is also at first order of smallness, so 
that r = er. As previously explained, the density and 
momentum is fully contained in f(°\ so that 
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Expanding the Boltzmann equation in this way and 
truncating the expansion to 0(e), we find 
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(12) 
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Gathering these terms according to their order of small- 
ness, we find 
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To derive the mass and momentum equations to the Eu- 
ler level, only the C(e°) terms are needed. The Navier- 
Stokes corrections to these equations are found by also 
including the 0(e 1 ) terms in the derivation. Similarly, 
it is possible to find further corrections at C(e 2 ) (known 
as the Burnett corrections) and beyond (super-Burnett), 
although these further corrections are negligible in prac- 
tical cases 0, Band have historically been viewed with 
some suspicion Q. 

Deriving the conservation equations from (|13p requires 
the second and third moments of /(°) [9]. In linearised 
form, these are 



i j ■ 



(14a) 



£,£,6./" < - Po(ui8 jk + u 3 5 lk + UkStj), (14b) 



where 5ij is the Kronecker delta. 

Taking the zeroth to second moments of (|13a|) . us- 
ing (|14p . and linearising, we find 



dp du t 



dt 



dx- 



du, dp 



(15a) 
(15b) 




duj ^ du k 
d Xi 13 dx k 



-S -- -U {1) 



(15c) 



where IT>V = J &£,7 (1) d£. With the O(e ) approxima- 
tion d/dto = d/dt, the two first equations are equivalent 
to linearised versions of the mass equation (|3|) and the 
Euler momentum equation with Sq and Si in the 
place of Q and Fi, respectively. Taking the zeroth and 
first moments of (|13a[) and linearising, we find 0(e) cor- 
rections to the above equations, 
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The sum (|15a[) + e (|16ap directly leads to the lin- 
earised mass equation, given below in (|20a[) . Similarly, 
()15b|) + e (|16b|) leads to the momentum equation, though 



the unknown tensor II 

(i) 



(i) 

ij 



Ujj can be explicitly related to the Navier-Stokes 



stress tensor using (|15c[) . The pressure time derivative 
can be rewritten assuming a nearly isentropic process and 
using (|15a[) . resulting in 

dp dp dp 2 / du k \ , , 

^ = ^ = H^%)' (17) 

Substituting for the speed of sound using ((9]), the diago- 
nal terms in (|15c[) become 



dp du k , 2 . du k 2 
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Thus, we find 
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+ T {Sij ~ CqSo) ■ 
The mass and momentum equations can now be ex- 
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plicitly found as described above, 



equation 



duj dp_ _ s _ p_ d(Sjj - SjjclSp) da-j 
dt dx t 1 p dxj dxj ' 

(20b) 

The momentum equation contains a deviatoric stress ten- 
sor 



dm diij 



da 



2 du k 



VB$ij^-, (21) 

dx k 



with shear viscosity fi = pqt and bulk viscosity \ib In = 
(5/3 — 7). This value of the bulk viscosity in the limit of 
rapid transfer of energy between translational and inner 
degrees of freedom has previously been found using more 
rigorous kinetic theory [l(| • 

Comparing (|20ap to the classical mass equation ([3]) , So 
appears in the place of the mass flux Q, which could be 
expected considering the interpretation of Sq as a mass 
flux. Comparing (|20b[) to the Euler-level momentum 
equation (0)), Si appears in the place of the body force, 
which had been neglected from the Boltzmann equation. 
Sij is also present inside a source term with a single spa- 
tial derivative and a small coefficient n/po in front. 



V. THE WAVE EQUATION 

The wave equation can be found as usual from the mass 
and momentum equations as <9 (|20ap /dt — <9 (|20b[) /<9z,. 
Using the isentropic relation (|17|) in the wave equation 
operator, and thus neglecting sound absorption due to 
thermal conduction and relaxation we find 
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(22) 



Comparing with (fT]) and © , we find a monopole strength 
To = dSo/dt, a dipole strength Ti = —Si, and a 
quadrupole strength which is not fully resolved but in- 
volves S^, S , and er^ . 

The deviatoric stress tensor o\' . can be resolved as 



d 2 u, 
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(23) 

The last term in the parenthesis contributes to sound 
absorption and is neglected in line with previous ap- 
proximations. The first parenthetical term contributes 
to the quadrupole strength. Using the property V 2 = 
Sij(d 2 jdxidxj), we find a fully resolved isentropic wave 
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H d 2 (Sij - 35ijp S /p ) 
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Finally, we find the quadrupole strength as Tjj — 
(fj,/p )(Sij - 3Sijp So/po). 

The coefficient fi/po in front of the quadrupole strength 
is typically on the order of 10 _10 s in gases. Its small 
magnitude means that the quadrupoles generated by the 
Boltzmann equation source term s tend to be negligible 
compared to the monopoles and dipoles. 



VI. SUMMARY AND CONCLUSION 

A source term s in the Boltzmann equation represents 
particles appearing or disappearing throughout the fluid 
with some distribution of particle velocities. As adding 
a source term to the mass equation allows modeling pul- 
sations of small bodies throughout the fluid [l[ , a source 
term in the Boltzmann equation would allow modeling 
more general vibrations of such small bodies. 

From this modified Boltzmann equation, the mass and 
momentum conservations equations (|20[) were derived un- 
der the common acoustic approximation of linearity and 
constant entropy. These equations gain source terms 
given by the moments (|10|) of s. The mass equation gains 
a source term Sq. To the Euler level, the momentum 
equation gains a source term Si, and to the Navier-Stokes 
level, it gains a source term involving Sij and Sij So- 

The wave equation derived from these two conserva- 
tion equations contains multipole source terms. The 
monopole strength is dSo/dt and the dipole strength 
is —Si- The quadrupole strength, which comes out of 
the Navier-Stokes level of the momentum equation, in- 
volves Sij and Sij Sq, but has a negligible magnitude. 
That the quadrupole strength is so much smaller than 
the monopole and dipole strength could be expected, as 
the Navier-Stokes level terms are one order higher in the 
small Knudsen number than the Euler-level terms. 

Similarly, going to the (9(Kn 2 ) Burnett level might 
lead to a d 2 Sijk/ 'dxjdxk term in the momentum equa- 
tion, leading to an octupole term in the wave equation. 
However, since this term would be at C(Kn 2 ), it would 
be even more negligible than the quadrupole term. 

Going back to the point of modeling general vibra- 
tions of small bodies throughout the fluid, this analysis 
indicates that the vibrations of such small bodies can ra- 
diate as monopoles and dipoles, but only very weakly as 
higher-order multipoles. 

Comparing this analysis with the analogous analysis 
for the lattice Boltzmann method [H, we find that the 
analogous source term in the lattice Boltzmann equation 
can only radiate quadrupoles effectively due to a fortu- 
itous discretisation error that occurs when discretising 



5 



the Boltzmann equation © in space and time using the zoidal method (H , which results in a scheme fully consis- 
first-order rectangle method. Discretising with the trape- tent with ([5]), would therefore also lead to a vanishing 

quadrupole term similarly to what has been shown here. 
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